Analytical results for long time behavior in anomalous diffusion 
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We investigate through a Generalized Langevin formalism the phenomenon of anomalous diffusion 
for asymptotic times, and we generalized the concept of the diffusion exponent. A method is 
proposed to obtain the diffusion coefficient analytically through the introduction of a time scaling 
factor A. We obtain as well an exact expression for A for all kinds of diffusion. Moreover, we show 
that A is a universal parameter determined by the diffusion exponent. The results are then compared 
with numerical calculations and very good agreement is observed. The method is general and may 
be applied to many types of stochastic problem. 
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I. INTRODUCTION. 

The study of systems with long range memory reveals 
some physical phenomena that are still not well under- 
stood, especially in systems which are outside the state 
of equilibrium or those in which the existence of anoma- 
lous diffusion is verified [T^S]. Here we show a simple 
analytical method which describes the behavior of the 
diffusion for large and intermediate times. In order to do 
that, we first generalize the concept of the diffusion ex- 
ponent. Then we present a conjecture to obtain, through 
the introduction of a time scaling factor A, an analyti- 
cal asymptotic result for the diffusion coefficient for long 
times. We obtain the scaling factor exactly and we show 
as well its universal behavior. We derive a numerical 
method to obtain the correlation function of velocities 
for an ensemble of particles from any given memory. We 
compare both methods and we obtain excellent agree- 
ment. The method has general application in the study 
of stochastic processes and it could be applied to several 
situations of physical interest. 



II. GENERALIZED LANGEVIN EQUATION 
AND DIFFUSION. 



The generalized Langevin equation (GLE) is a stochas- 
tic differential equation which can be used to model sys- 
tems driven by colored random forces. For the velocity 
operator v(t) this equation can be written as, 



m / T(t - t')v(t')dt' +f(t), (1) 



dv(t) 
dt 



where T(t) is the retarded friction kernel of the system, 
or the memory function. Here, is a stochastic noise 
subject to the conditions = 0, (£(t)v (0)) = 0, and 



where C^(t) is the correlation function for and the 
angular brackets denote an average over the ensemble of 



particles. Equation (12 
theorem (FDT) $\ M 



is Kubo's fluctuation dissipation 
. The presence of the kernel T(t) 
allows us to study a large number of correlated processes. 
In the real world, the vast majority of problems are non- 
Markovian, i.e., there is correlation between the various 
stages of dynamic evolution. This property is what we 
call memory, and it makes remote events of the past im- 
portant to dynamic events in the present time. 

Using the GLE it is possible to study the asymptotic 
behavior of the second moment of the particle movement, 



lim (x'(t)} 

t—^OO 



2D(t)t ~ t c 



(3) 



which characterizes the type of diffusion presented by 
the system. Here, D(t) is the diffusion coefficient as a 
function of time. Moreover, for an asymptotic behavior 
of the form 



lim (x 2 (t)) ~ t a [\n(t)} 



±i 



(4) 



QWH£W£(o)) = mV«)r(t) 



(2) 



we shall define respectively an a diffusive behavior. 
Here the exponent a — a arises in analogy with the 
critical exponents in a phase transition. For example, 
in the two-dimensional Ising model the critical exponent 
for the specific heat is a = + because it does not have 
a power law behavior; rather it has In \T — T c \ behavior 
for temperatures T close to the transition temperature 
T c . This generalized nomenclature is pertinent here since 
there are quite a large number of possibilities of combi- 
nations for logarithmic and power law behaviors. 

In this way the behavior of D(t) can be determined 
using 



lim D(t) = lim lim / C v (t') exp(-zt')dt' = limR(z), 

t— >oo i— >oo z— >0 Jq z— »0 

where R(t) = C v (t)/C v (0), with C„(0) = 1, and R(z) is 
the Laplace transform of R(t). For t — > oo and normal 



2 



diffusion this is the Kubo formula [TO]. The limits can 
be justified using the final value theorem (FVT) for a 
Laplace transform i.e., for any function g(t) with 
Laplace transform ~g(z) then lim^oo g(t) = lin^^o zg[z). 
Now a Laplace transform of the integral gives D(z) = 
R{z)/z, and we end up with the equation above. 

Now we multiply Eq. (IT]) by v(0) and take the average 
over the ensemble, with (f (t)w(O)) = 0, to obtain a self- 
consistent equation for R(t) in the form 



R(t) 



F(t - t')R(t')dt'. 



We then Laplace transform Eq. ^ to get 

1 



R(z) = 



z + T(z) 



(6) 



(7) 



Time correlation functions play a central role in non- 
equilibrium statistical mechanics in many areas, such 
as the dynamics of polymeric chains [T2], metallic liq- 
uids [13], Lennard- Jones liquids [2], ratchet devices [TS] 
IT5] . spin waves [T7], Heisenberg ferromagnets and dense 
fluids [18 . Consequently to invert this transform, or a 
similar one, is crucial. Unfortunately, in most cases it is 
not an easy task. In those situations, the use of numerical 
methods is an alternative to overcome this problem. Our 
main objective here is to show a process to obtain the 
asymptotic behavior analytically. Although the method 
can be applied to several situations, we concentrate here 
on the analysis of diffusion. 



III. THE CONJECTURE 

We claim that after a "transient time" r, i.e., for t > r, 
the leading term for D(t) will fulfill Eq. ^ within a given 
approximation. In this context t — > oo is equivalent to 
t ^> t. Now we imposed the scaling 

z -> X/t. (8) 

In order to determine A we rewrite Eq. ([5]) as 



lim D(t) = lim R(z = X/t) = lim -— , (9) 

t— ¥OD t—¥OD t—^-OO fit) 



where 



f(t) = x + tr(x/t). 

The derivative of Eq. ([9| yields 



(10) 



lim R 1 (t) = lim ±D(t) = lim [1 - 4 In [/(*)]]//(*)> 

t-+oo t-too dt t-voo clt 

(11) 

(12) 



while from the FVT we have 



lim J?2(i) = hm zR(z) = lim -77-r. 



The relative diference 

AR[t) = ^rV 1 = [x ~ l + t Jt ln [fm,x (13) 

should evolve to zero as t — > oo. For A ^ this yields 
the exact value 



A = l- lim t- In [/(«)]. 
t-yoo dt 



(14) 



The scaling works as long as the GLE, Eq. J7j), works. 
To obtain A we need more information about T(z), which 
may be different for every system. However, since our 
interest is in the asymptotic behavior, we can expand 
r(z) , in Taylor or Laurent series around z = 0, in the 
form 



T(z) - z"[a- 61n(z) - c/ ln(z)], 



(15) 



where a, 6, and c are positive constants. Note that we 
give especial attention to ln(z), since it will give us the 
behavior pointed out in Eq. Q. For 6 = this gives a 
diffusion with exponent a; for i) ^ this gives an a - , 
and for a = b = and c ^ we get an a + diffusion. 
If r(z) has another contribution, besides ln(z), that can 
not be expanded at the origin we keep it and expand 
the other parts. However, most of the memories in the 
literature can be cast in the form Eq. (151 for small z. 
Now we introduce Eq. (151 into Eq. ( 14 ) to obtain A = v 
for v < 1, and A = 1 for v > 1. Notice that it does not 
depend on a, 6, or c, which suggests a universal behavior. 

In our conjecture some points deserve attention: First, 
we are considering integrals, of the form Eq. ([5|, where 
the function R(t) is well behaved, and limited to — 1 < 
R(t) < 1, since C v (t) < C„(0). R(t) is such that it always 
has a well-defined behavior for finite t, even when the in- 
tegral diverges as t — > oo, as in superdifusion. Second, 
D(t) must have a leading term as t — ¥ oo, which deter- 
mines the diffusion. For example, the inverse Laplace 
transform of R(z) is 



R(t) 



1 

2tt7 



-\-ioc-\-rj 



R(z) exp(zt)dz. 



(16) 



— ioo+77 



Here the real number 77 is such that all the singularities lie 
at the left of the line joining the limits. Consider now Eq. 
( 15 1 with b = c = 0, and v < 1; then lini^o R(z) ~ z~ v , 



and 



/+ioo+r/ 
s~ u exp(s)ds cx t u ~ l , (17) 



where we have done the transformations s — zt and rj' = 
r\jt. For v > the only pole is at s — , and the 
condition in ?/ will be automatically satisfied. Now by 
direct integration on Eq. ^ we obtain D(t) cx t v . From 
the scaling we get the equivalent result 



lim D(t) = lim R(z = X/t) 



lim R(X/t) - t v . 



(18) 
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Note that the above exact result is not only for power 
laws, but for any function behaving as a power law for 
large t. We confirm as well the relation a = v + 1, ob- 
tained by Morgado et al. [1 j. Our results can be readily 
expressed as 



A 



f 



1 



—1<v<1, 

V > 1. 



(19) 



The factor A depends only on the diffusion exponent a, 
consequently it is universal. Moreover it will be the same 
for a or a*. For normal diffusion a = 1, or for a = l^, 
A = 0. However, we still can obtain the final value. Con- 
sider as example the Langevin equation without memory; 
for that we have R(t) = exp (—yt) and R(z) = (7+ z) -1 . 
From Eq. ([9]) we get 



lim D(t) = lim i?(A/t) = 



t 



'jt + X 



= 7 



while direct integration gives 



lim 



D(t) = lim 

t— >oo 



R(t')dt' = 1 - 



(20) 



(21) 



In this case the scaling yields correctly the wanted final 
value. 

Equation ([6]) imposes as well some requirements on 
R(t). First its derivative must be null at the origin, i.e., 
the integral in the right hand side must be null at t = 0. 
This is true except for non analytical memories, such as 8 
functions. Indeed, we do not expect exponential behav- 
ior of the form R(t) = exp (— -y | £ | ) with a discontinuous 
derivative at the origin • Second, in Eq. (I), for 

a bath of harmonic oscillators the noise can be obtained 
as [12] 

f(f) = / ^2k B Tg(oj) cos[ujt + (f>(u;)]du, (22) 



where < </>(w) < 2tt are random phases and g(oS) is the 
noise spectral density. The FDT yields 



r(t) 



<?(w) cos(ujt)du. 



(23) 



This shows that the memory is an even function of t. An 
analytical extension of T(z) in the whole complex plane 
has the property T(—z) = — T(z). Consequently, from 
Eq. {7}, R{-z) = -R(z), or R{-t) = R(t). In short, 
it requires well-behaved functions and derivatives. Even 
functions have zero derivatives at the origin as required 
before. 



IV. THE BALLISTIC DIFFUSION 

Let us consider the spectral density 
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FIG. 1. Correlation function R(t) as a function of time t. 
We use the memory ( |25[ ) and numerical integration. Curve a, 
(jj s = 1, and 6 = 1; curve b, u 3 = 5, and b = 1/2. 



This is a generalization of the Debye density of states. 
Here b > is a dimensionless constant, and oj s is a cutoff 
frequency. For j3 ^ we get anomalous diffusion. In 
particular, for j3 = 1 we introduce Eq. (24) into Eq. (23) 
to obtain 



2 / sin(w s t) 



r(t) = buj 

with the Laplace transform 



cos(w s £) — 1 



I» = yln 



(25) 



(26) 



First, we have the analytical function D(t) = R(z = A/t); 
second, from Eq. (14) we obtain lim^oo A = 1, exactly. 
This is ballistic diffusion of the form a — 2~. 

Now we compare the analytical asymptotic with a nu- 
merical solution of Eq. (J6|. To do this, we rewrite this 
equation in a discrete form, and then we expand it up to 
terms of order At 2 " to obtain 



R(t + At) = R(t -At) + 2J2 R (2k ~ l) {t) 



(At) 



2fc-l 



fe=0 



(2k- 1)! 



(27) 



where R( n \t), is the time derivative of R(t) of order n. 
Note that this expansion eliminates all the even deriva- 
tives. Now we can obtain all R(t + At) from the sequence 
of the previous value of R(t), starting from R(0) = 1. 
From these values, its possible to get the diffusion coeffi- 
cient through direct integration of Eq. ^ . 

In Fig. [l] we plot the correlation function R(t) as a 
function of time t. The curves correspond to the numer- 
ical solution, and are calculated using Eq. (27), and Eq. 
(25 ) with At = 10 -5 . For curve a, u> s = I, and 6 = 1; for 



curve 6, oj s = 5 and 6 = 1/2. 

In Fig. [2] we plot the diffusion coefficient D(t) as a 
function of time t. The oscillatory curves corresponds to 
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[(l + a)z] _1 or R(t) = [l + a] -1 , and by direct integration 
we get D(t) = t/(l + a) exactly. This is ballistic a = 2 
diffusion. If we apply Eq. ([9| we obtain the same result 
with A = 1. Since from the relations (15) and (19) the 



FIG. 2. Diffusion coefficient D(t) as a function of time t. 
Curve a, uu a — 1, and 6=1; curve b, u s = 5, and 6 = 1/2. 
The oscillatory curves are the numerical result. The curves 
without oscillations are the analytical asymptotic limit. We 
see in curve b that the two curves collapse onto a single one. 



the numerical solution and are calculated from the data 
of Fig. 1. The curves without oscillations correspond 
to the analytical asymptotic limit, Eq. ([9|, with memory 
Eq. ( 26 1 . Here we see that the asymptotic curves are 



mean values of the oscillatory ones. In this range the fit 
yields For curve a, A = 0.928 ± 0.002, and for curve b, 
A = 0.948 22 ± 0.000 01. We see in curve b that the two 
curves collapse onto a single one. Here the transient time 
r to which we refer before Eq. ^ is a decreasing function 
of b/oj s . The value of A approaches the exact value 1 as 
the ratio b/uj s decreases, or as time increases. This shows 
the efficiency of the scaling; even before convergence is 
fully established, curve a, the asymptotic curve gives us 
an average value that can be used to understand the main 
characteristics of the process. 

Consider now T(z) = az, exactly. That means R(z) — 



value of A does not depend on In (z), this result is exactly 
what we get from Eq. ( |26| . There are important differ- 
ences between the a = 2~ diffusion, which according to 
the Khinchin theorem [51 [2T] is ergodic, and the a — 2 
diffusion, which does violate ergodicity. This distinction 
was not possible before the generalization of the diffusion 
exponent we present here. 



V. CONCLUSION 

In this work we generalize the concept of the diffusion 
exponent, and we propose a conjecture to investigate the 
asymptotic limits of anomalous diffusion, through the in- 
troduction of a time scaling factor A. We obtain the scal- 
ing parameter exactly and we show that it is universal 
and depends only on the diffusion exponent. We analyze 
the ballistic diffusions a — 2~ and a = 2, both ana- 
lytically and numerically. The method can be useful as 
well to analyze large amounts of data in stochastic pro- 
cesses [5] , and in different fields of science where is nec- 
essary to inverse a Laplace transform of the form of Eq. 
([7]). The phenomenon of diffusion also poses challenges 
in the understanding of fundamental concepts in statis- 
tical physics, such general properties as the correlation 
function [19] , ergodicity j6[ HI l2~iTf23] , and the Khinchin 
theorem 6, 21 . In nonlinear phenomena, such as growth 
and etching [24] . analytical results are rather difficult to 
obtain. In this way we hope that this work may inspire 
research into similar asymptotic limits. 
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